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On the response of quantum linear systems to 
single photon input fields 

Guofeng Zhang Matthew R. James 



Abstract 

The purpose of this paper is to extend linear systems and signals theory to include single photon 
quantum signals. We provide detailed results describing how quantum linear systems respond to multi- 
channel single photon quantum signals. In particular, we characterize the class of states (which we call 
photon-Gaussian states) that result when multichannel photons are input to a quantum linear system. 
We show that this class of quantum states is preserved by quantum linear systems. Multichannel photon- 
q-( Gaussian states are defined via the action of certain creation and annihilation operators on Gaussian 

states. Our results show how the output states are determined from the input states through a pair 
of transfer function relations. We also provide equations from which output signal intensities can be 
■ computed. Examples from quantum optics are provided to illustrate the results. 

Index Terms — Quantum linear systems; continuous mode single photon states; Gaussian states 
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I. Introduction 

One of the most basic aspects of systems and control theory is the study of how systems 
respond to input signals. Well known tools, including transfer functions and impulse response 
functions, allow engineers to determine the output signal produced by a linear system in response 
to a given input signal. Such knowledge is needed for engineers to enable them to analyze and 
design control systems. 

As is well known, signals other than deterministic signals are important to a wide range 
of applications. There is a well-developed theory of dynamical systems and non-deterministic 
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signals, within which Gaussian signals play an important role. Indeed, linear systems (with Gaus- 
sian states) driven by Gaussian input signals provide the foundations for Kalman filtering and 
linear quadratic Gaussian (LQG) control, as well as many other developments and applications, 
ED, (22), [0, fl2J, [EU- It is well known that for a classical linear system G initialized in a 
Gaussian state and driven by Gaussian white noise process w(t) satisfying E[u>(t)] = and 
K[w(t)w(t') T ] = S(t - t'), 

x(t) = Ax(t) + Bw(t), x(t ) = x, (1) 

y(t) = Cx{t)+Dw(t), 

where E denote mathematical expectation, the mean x(t) = ~E[x(t)] and covariance E(t) = 
E [(x(t) — x(t))(x(t) — x(t)) T ~\ satisfy the differential equations 

x(t) = Ax(t), (2) 
£(t) = AT J (t) + T,(t)A T + BB T . (3) 

These equations characterize the dynamical evolution of the Gaussian distributions of the state 
variables x(t). Expected values of quadratic forms can easily be evaluated, for instance in the 
zero mean case F,[x T (t)Mx(t)] = Tr[M T S(t)]. In the frequency domain, the spectral density 
R out [iu] of the output process y(t) is related to the spectral density R in [iuS\ of the input process 
w(t) via the transfer relation 

R out [iu] = E G [iu]R in [iu]E G [iu]^ , (4) 

where Eg is the transfer function for the system G. Differential equations of the types © and 
®, and transfer relations like ©, play fundamental roles in classical linear systems and signals 
theory, ffl, fl22l, flU, 0, El. 

Quantum linear systems are a class of open quantum systems fundamental to quantum optics 
and quantum technology, 0, J20l, [HU, El, 03, (301. The equations describing quantum 
linear systems (see Section III-AI below) look formally like the classical equations (OQ), but they 
are not classical equations, and in fact give the Heisenberg dynamics of a system of coupled open 
quantum harmonic oscillators, [|43l . |fT9l . These quantum systems are driven by quantum fields 
that describe the influence of the external environment (e.g. light beams) on the oscillators. 
In quantum optics the fields play the role of quantum signals. As pointed out in the paper 
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[301, a large fraction of quantum optics literature concerns fields in coherent states (a type of 
Gaussian state), and as a consequence quantum optical systems driven by coherent fields are 
well understood. Indeed, when a quantum linear system, initialized in a Gaussian state, is driven 
by a Gaussian field, the state of the system is Gaussian, with mean and covariance satisfying 
equations of the form © and ©■ This fact is of basic importance to Gaussian quantum systems 
and signals theory, and for example has been exploited for the purpose of quantum H°° and 
LQG control design, Ea.Ba.lMED.Ea.ltt 

In recent years, due to their highly non-classical properties, single-photon light fields have 
found important applications in quantum communication, quantum computation, quantum cryp- 
tography, and quantum metrology, etc., [|20l . 0, OOl . Unlike Gaussian states, a light pulse in 
single-photon state contains one and only one photon, and is thus highly non-classical. While 
much of the optical quantum information literature deals with single modes (namely discrete 
variables like polarization) of light and static devices (like beamsplitters), the importance of 
continuous mode photons and dynamical devices is becoming clear, [|35l . [|29l , ll30l . However, 
single photon states of light are not Gaussian, and so the relatively well developed quantum 
Gaussian systems and signals theory is not directly useful for quantum optical systems driven 
by single photon fields. 

The purpose of this paper is to extend linear systems and signals theory to include single 
photon quantum signals. We build on the results in [|29l to describe how quantum linear dynamical 
systems respond to multichannel continuous mode photon fields from a system-theoretic point of 
view. We show, for example, that when a quantum linear system G with no scattering (equations 
(fT4l)-([T5T) below with S = I) is driven by multichannel photon fields, the mean a(t) = Tr[pd(t)] 
and covariance E(t) = Tr[p(a(t) — a(t))(a(t) — a(t))t], where p is the initial joint system-field 
density operator, satisfy the differential equations 

&(t) = Al(t), (5) 
t(t) = AE(t) + E(t)A* + Brt(£(*)) + r(£(t))Bi + BFBl (6) 

In equation ©, F is a matrix depending on the Ito products of the input fields ll32ll . [fT3l . and 
r(£(i)) is a matrix depending on the pulse shape matrix £(£). Equation © is crucial to the study 
of intensity of output fields, cf. Sec. ITTll 

When multichannel photons are input to a quantum linear system, the output state can be quite 
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Fig. 1. A quantum linear stochastic system G driven by multichannel field in a photon-Gaussian state (Definition [TJ. The 
multichannel input and output field states are denoted Pi in ,R in and P( out ,R out respectively. The system G transfers the input 
pulse shape £j n and covariance function Ri n to the output pulse shape £ ou t and covariance function Rout respectively. In this 
paper it is assumed that the system is initially isolated from the input field. 

complex. To accommodate the types of states that can be produced from multichannel photon 
inputs, we define a class T of quantum states, which we call photon-Gaussian states. A state 
(density operator) p G T is specified by a matrix £(t) of functions (or pulses, a multichannel 
generalization of wavepackets), and a Gaussian spectral density R[iu). We sometimes express 
these states as p^R, as in Fig.Q] Our main result, Theorem [51 states that if the photon-Gaussian 
state P£ in ,n in G T is input to a quantum linear system G initialized in the vacuum state, then the 
steady-state output state is also a photon-Gaussian state P( ouU R out G T . Moreover, the transfer 
Ri n H- Rout is given by the above relation ©, while the transfer £ in H- £ out is given by 



where E G is the transfer function of the quantum linear system G defined in Sec. [Til This result 
provides a natural generalization of the well-known Gaussian transfer properties of classical linear 
systems to an important class of highly non-classical quantum states that includes single photon 
states. Results of this type are anticipated to be of fundamental importance to the analysis and 
design of quantum systems for the processing of highly non-classical quantum states. While 
most of this paper is concerned with questions of analysis, we include a short section on 
synthesis, generalizing the work ll29ll . Here, a quantum linear system is designed to manipulate 
the wavepacket shape of a single photon. This is an example of coherent control, ||29l . 

Notation. Given a column vector of complex numbers or operators x = [ x\ ■ ■ ■ x k } T where 
k is a positive integer, define x # = [ x\ ■ ■ ■ xl ] T > where the asterisk * indicates complex 
conjugation or Hilbert space adjoint. Denote x^ = (x*) T . Furthermore, define the doubled-up 
column vector to be x = [ x T (x # ) T ] T - The matrix case can be defined analogously. Let I k be 
an identity matrix and 0^ a zero square matrix, both of dimension k. Define = diag(ifc, —Ik) 



£ out [s] = Z G [s]£ in [s}, 



(7) 
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Fig. 2. A quantum linear stochastic system G with input b and output b out 



and 0fc = [ Ik] —Ik ] (The subscript "k" is omitted when it causes no confusion.) Then 
for a matrix X 6 C 2:jx2k , define X b = J^X^Jj. ® c denotes the Kronecker product, m is the 
number of input channels, and n is the number of degrees of freedom of a given quantum 
linear stochastic system (that is, the number of oscillators). \<p) denotes the initial state of the 
system which is always assumed to be vacuum, |0) denotes the vacuum state of free fields. 
Given a function f(t) in the time domain, define its two-sided Laplace transform [23, Chapter 
10] to be F[s] = Sfb{f(t)}(s) := e~ st f(t)dt. When s = ico, we have the Fourier transform 



FUlo} 



i: 



-iuit 



f(t)dt. Given two constant matrices U, V G C rxk , a doubled-up matrix 



A (U, V) is defined as 



A(U,V) 



U V 

V* u* 



(8) 



Similarly, given time-domain matrix functions E (t) and E + (t) of compatible dimensions, define 
a doubled-up matrix function 

E~(t) E+(t) 



A(E~(t),E + (t)) := 
Then its two-sided Laplace transform is 

A(E-[s},E+[s}) = J? b {A(E-(t),E + (t))}(s 



E+(t)# E~{t)* 



(9) 



(10) 



E-[s] E + [s] 
E + [s*}* E~[s*}* 

Finally, given two operators A and B, their commutator is defined to be [A, B] = AB — BA. 



II. Quantum Linear Systems 

In this section quantum signals and systems of interest are introduced, Fig.[2l Quantum systems 
behave in accordance with the laws of quantum mechanics, and in Sec. Ill- Al we summarize the 
dynamics of a quantum linear system driven by external quantum fields. These models feature 
inputs and outputs, corresponding, for example, to light incident on and reflected by the system, 
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|fT3l , 11431 . In Sec. III-BI we write down explicitly the input-output relations, using a notation for 
the impulse response motivated by physical annihilation and creation processes. Since one of 
our objectives is to study the steady-state response of the quantum linear system, we present in 
Sec. III-CI the steady-state versions of the input-output relations, as well as the transfer function 
defined in terms of the two-sided Laplace transform. Our analysis will require the stable inversion 
of the transfer function, and this is presented in Sec. III-DI Sections III-EL HI-FI and III-GI provide 
the definitions of Gaussian and photon states needed in this paper. Sec. HI-HI then describes how 
output means and covariances are determined from the corresponding input quantities. 

A. Dynamics 

The open quantum linear system G, shown in Fig. [2l is a collection of n interacting quantum 
harmonic oscillators a = [a 1; . . . , a n ] T (defined on a Hilbert space Sj G ) coupled to m boson fields 
b{t) = [61 (t), . . . , b m {t)] T (defined on a Fock space g) [El, flU, 02. Here, aj (J = 1, . . . , n) 
is the annihilation operator of the jth quantum harmonic oscillator satisfying the canonical 
commutation relations [aj,^] = 5j k . The vector b(t) represents an m-channel electromagnetic 
field in free space, which satisfies the singular commutation relations 

[&,■(*),&*(*')] = W(*-0, [&i(*)A(01 = [&,-(*), &£(*')] = 0, j,k = l,...,m, VM' e R. (11) 

The operator bj(t) (j — 1, . . . , m) may be regarded as a quantum stochastic process; in the case 
where the field is in the vacuum state (denoted |0), 11321 . tfT3l . Il43l ). this process is quantum white 
noise. The integrated field operators are given by B(t) = f*b(r)dr and B*{t) = f* b#(r)dr, 
which are quantum Wiener processes. The gauge process is given by A(t) = J b*(r)b T (r)dr = 
(Aj k (t)) ■ k=1 m with operator entries A jfe on the Fock space T . Finally in this paper it is assumed 
that these quantum stochastic processes are canonical, that is, they have the following non-zero 
Ito products 

dB^dBKt) = 5 3k dt, dA jk dB!(t) = S kl dB*{t), (12) 
dBjtydkuty) = 6 jk dBi(t), dk jk (t)dk lm (t) = 6 kl dA jm (t), (j, k, I = 1, . . . , m). 

The system G can be parameterized by a triple (S-,L,H). In this triple, S- is a scattering 
matrix (satisfying SLSl = SLS- = I m ). The vector operator L is defined as L = C^a + C + a#, 
where C_ and C + £ C mXn . The Hamiltonian H = |a^A f2 + ) a describes the initial internal 
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energy of the oscillators, where fi_,f2+ G C nxn satisfy f2_ = fit and Q + = fiT. With these 
parameters, the Schrodinger's equation for the system G (with initial internal energy H) and 
boson field is, in Ito form ( |fT3l Chapter 11]), 

dU{t, t ) = ^Tr[(S- - I m )dA T } + dB\t)L - L ] S^dB(t) - (^L j L + iH)dtj U{t, t ), t> t , 

(13) 

with U (to, t ) = I (the identity operator) for all t < t . Moreover, by means of Ito rules, it can be 
shown that the operator U(t, t ) satisfies d(U(t,t )U*(t,t )) = d(U*(t,t )U(t,t )) = 0, Vt > t , 
and since U(to,t ) = I we see that U(t,to) is unitary. 

In the Heisenberg picture, the system operators evolve according to d(t) = U(t,to)*dU(t,t ) 
(componentwise on the components of a). The output field b ou t{t) (which carries away infor- 
mation about the system after interaction) is defined by b ou t(t) = U(t,to)*b(t)U(t,to) (compo- 
nentwise on the components of b(t)). Consequently, by Eq. (fT3l) . the dynamical model for the 
system G can be written as 

a(t) = Ad{t) + BSb (t) , d(t ) = a, (14) 
kut(t) = Cd(t) + Sb(t), (15) 

in which system matrices are given in terms of physical parameters by S = A (S^O) , C = 
A (C_, C + ) , B = -C\ A = -\&C - iJ n A Q + ). 

Remark 1: Eqs. (fl4)) and (TTST) are quantum linear systems, which can be obtained using the 
definitions a(t) = U*(t,t )a(to)U(t,t ) and b out (t) = U*(t,t )b(t)U(t,t ) and making use of 
Ito rules and the commutation relations for annihilation and creation operators, ll43l . 

The system G is said to be asymptotically stable (equivalently, exponentially stable) if the 
matrix A is Hurwitz H4l Sec. Ill- A]. More details on quantum linear stochastic systems can be 
found in, e.g., Q31, 03, EH, MB and references therein. 

Remark 2: Quantum linear systems have been widely used in quantum optics. Moreover, 
They have also been used in opto-mechanical systems, e.g., Eqs. (15)-(18) in the Supplementary 
Information of Reference ll27l and Eqs. (9)-(10) and the line below Eq. (10) in Reference ll39ll . 
They also appear in circuit quantum electrodynamics (circuit QED) systems, e.g., Eqs. (18)-(21) 
in Reference [28 j. 
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B. Input-output relations 

The output may be expressed in terms of the input and initial system variables 

b out (t) = Ce A{t - to) a + f Ce A(t ~ r) BSb{r)dr + Sb(t), t > t . (16) 

We find it convenient to express this input-output relation in terms of impulse response functions. 
Define 



9G-{t) 



6(t)S. 
0. 



-c. 



t > 
t < 



■[ C_ C + }e 



At 



-cl 



s 



# 



t > 



t < 



(17) 



9a(t) 



(Later we use g J Q_ and (f G+ (j,k = 1, . . . , m) to denote the entries of g G - and g G + on the jth 
row and /cth column, respectively.) The impulse response function for the system G is 

5{t)S- Ce At C^S, t > 
0, t < 

It can be checked that ga(t) defined in Eq. (fT8l) is in the form of 

g G (t) = A(g G -(t),g G+ (t)). 

Therefore the input-output relation (fT6l) may be expressed in a more compact form 

kutit) = Ce^-^a + f g G {t - r)b(r)dr. 



(18) 



(19) 



(20) 



to 



C. Steady-state input-output relations 

Assume that the system (fl4l) is asymptotically stable. Letting t — > — oo and noticing ([TBI . 
Eq. (f20l) becomes 



b ut(t) 



(21) 



5f G (i - r)b(r)dr. 

Let Hg[s], S g -[s], S g +[s], b[s] and feoutfs] denote the two-sided Laplace transforms of g G (t), 
9G-{t), 9G+{t), b(t) and b out {t), respectively. Then by the above definitions and standard proper- 
ties of the two-sided Laplace transform we have the transfer function relation b out {s] = E G [s]b[s], 
where 

S G [s] = A(S -[s],S 0+ [s]). (22) 
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D. Flat-unitary property and stable inversion 

For later use we record here some important inversion results for quantum linear systems. 
Proposition 1: The reciprocal of the transfer function S G [s] is given by 

Zeis]" 1 = ~ G [- S f = A(S G -[- S *]t,-H G+ [- S ] T ). (23) 

Consequently, we have the fundamental flat-unitary relation 

~ G H b ~ G H = ~ G M2 G M b = hm, Vu6l. (24) 

Proof: Equation d23l follows from lfT6l Eq. (74)], and the flat-unitary relation (f24|) follows 
on setting s = iui. ■ 
We now apply the results and methods of ll38ll to find a stable inverse of the quantum linear 
system G. Define 

g a -x(t) := ^{Eais]- 1 }®, (25) 

where Jz^ 1 is the inverse two-sided Laplace transform If23l Chapter 10]. 
On the basis of Proposition \T\ the following result can be established. 

Lemma 1: Assume that the system G (Eq. (fT4"l) ) is asymptotically stable. Then the impulse 
response of the stable inverse of G is given by 

g G -,{t) = A {g G -{-t)\ -g G+ (-t) T ) = I °\ . * > ° ■ (26) 

Proof: By definition of the two-sided Laplace transform, we have Jfb{f(—t)*}(s) = f[—s*]* 
and <£?;,{/(— £)}(s) = /[— s]. Combining this with expression (|23T) in Proposition [Q we obtain 
the first equality in d26l) . The second equality in (l26l) follows from the definitions of the impulse 
responses ga-(t) and g G +(t) in terms of the system matrices. ■ 

E. Gaussian system (single mode) states 

A stationary Gaussian system state p s on f) G is specified by its characteristic function lfl4l . 

Tr[p s exp(^ f a)] = exp(-^ t S^ + iz ] (3), Wz e C n , (27) 

where f3 = a(t) = Tr[p s a(t)], and E = Tr[p s (a(t) — f3)(a(t) — /3) t ] is a non-negative Hermitian 
matrix. In general, S has the form 

r In + N T M 

(28) 

Aft N 
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In particular, the ground or vacuum state \<j>) is specified by (5 = and E 9 



In 
0„ 



F. Gaussian field states 

Depending on the nature of the boson fields, input signals b(t) can be in various states. In this 
section we consider m-channel Gaussian field states p g , |[T4|. Given a function / G L 2 (M, C m ), 
define an integral functional £H on the Hilbert space L 2 (R, C 2m ) to be 



W)(t) == 



R g (t,r)f(r)dr, 



(29) 



where R g (t,r) = R g (r,ty > 0. The m-channel Gaussian field state can be specified by the 
characteristic function 



r»OC 

Tr[p 9 exp(z / f{tfb{t)dt)) = exp(- -(/>/) + i(m\f)), 



(30) 



where m e L 2 (R, C m ), and (-|-) denotes the inner product in the Hilbert space L 2 (M, C 2m ). It 
can be checked that the mean value is Tr[p g b(t)] = m(t) and the covariance function is 

Ti[p g (b(t)-m(t))(b\r)-xh(ry)] = R g (t,r), Vr,t G R. (31) 

In general, the covariance function R g (t,r) has the form 

S(t - r) + u(r, t) T p{t,r) 
p(r,tf u(t,r) 

for real r and t. When p g is stationary, R g (t,r) depends on the difference between t and r, 
instead of their particular vales. In this case, we may use R 9 (t) to replace R g (t, r). In particular, 
the quantum vacuum field state |0) is specified by 



R g {t,r) 



(32) 



Ro(r) = S(t) 



I m 
m 



(33) 



G. Photon field states 

Now we introduce another type of field states: the continuous-mode single photon pure states. 
In the one channel case we denote it by |\&) = |1„), as defined in EH Eq. (6.3.4)], Il2~9l Eq. (9)] 

\l v )=B*(v)\0), (34) 
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R(t,r)=E lv [b(t)#{r)] = 6(t-r) 



1 













+ 



(35) 



where |0) is the vacuum state of the field as defined in Sec. III-FL and B*{v) := u{r)b*{r)dr. 
Here, v is a complex-valued function such that \u{r)\ 2 dr = 1. Expression (l34l) says that 
the single photon wavepacket is created from the vacuum using the field operator B*(y). Using 
the relation b(t)\\ v ) = u(t)\0), we see that the operator B(u) := J™ u(r)*b(r)dr annihilates a 
photon, resulting in the vacuum: B(y)\l v ) = |0). 

The field operator b(t) is a zero mean quantum stochastic process with respect to the single 
photon state. The covariance function is given by 

u(r)*v(t) 
v{r)v{t)* 

The gauge process A(t) (recall Section III-AI) for a single photon channel takes the form 
A(t) = J^n(r)dr, where n(t) = b*(t)b(t) is the number operator for the field. The intensity 
of the field is the mean n{t) = (l^|rz(t) 1 1^) = |^(t)| 2 , an important physical quantity that 
determines the probability of photodetection per unit time. 

H. Mean and covariance transfer 

In this section some basic covariance transfer results for the quantum linear system G are 
presented, which are the quantum counterparts of the well-known classical results, e.g. [22, Sec. 

I. 10.3], but adapted to take into account the non-commuting system and field variables. 
Consider a quantum linear system G initialized at time to in a state p = \(f>){(f>\ ® p field, where 

|0) is the vacuum system state and the field state pjieid satisfies Ti[pf ie idb{t)] = (particular 
choices of Pfi e id will be made below). By taking expectations in Eq. (fl4l) . we find that the 
mean a(t) = Tr [pa (£)] satisfies Eq. §5$ since the field has mean zero. Also, since the system is 
initialized in the ground state |</>), we have d(t) = a (to) = for all t > t . Define the matrix 

I\t) = Tr[pa(tp(t)}, (36) 

and a matrix F by 

Fdt = Tr[p field dB(t)dB\t)]. (37) 

Lemma 2: The quantum linear system G initialized in the state p (ground system state and 
zero mean field state) has zero mean a(t) =0 for all t > t and the covariance matrix £(£) = 
Tr [pa (t) at (t)] satisfies the differential equation 

±{t) = AE{t) + E{t)Ai + BST\t) + T{t)S j B^ + BSFS j B\ (38) 
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with initial condition £(£ ) = where F is given in Eq. (1371) . 

Proof: The proof of Eq. (l38l) follows by taking expectations of the differential 



d{a{t)a\t)) = (da(t))a\t) + a{t)da\t) + da(t)dtf{t) (39) 
= Ad{t)d\t)dt + BSdB(t)d\t) + d{t)d\t)A^dt + d{t)dB\t)S ] B ] + BSFS ] B ] dt 

and noticing that Ti-[pdB(t)at(i)] = Tr[pa(t)dBt(t)]t. ■ 

Note that expected values of quadratic forms may easily be evaluated in terms of £(t), for 
example, Tr\pa){t)Ma{t)] = Tr[Af T S(t) T ] -Tr[Af r J]. 

Remark 3: Further information regarding the dynamics of T(t) will be given in Section ITlI-B I 
for the case of multichannel photon fields. 

Suppose now b out {t) is the steady-state output field defined by (121) . Define the input and 
output covariances R in (t,r) : = Tr[p6(t)6 t (r)] and R mt (t,r) := Ti\pb out (t)b ] out {r)}. 

Theorem 1: Assume that the system (fT4l) is asymptotically stable. Let the input field have 
covariance Ri n (t,r). Then the steady-state output covariance is given by 



/OO rOO 
/ 9c(t - r x )R in (T 1 , T 2 )g G (r - T 2 ) ] dTidT 2 . 
-OO J —OO 



(40) 



Now suppose that b(t) is stationary with respect to the field state pfieid- Write R[iu>] for the 
spectral density matrix (where R[s] is the two-sided Laplace transform of R(t)). 

Theorem 2: Assume that the system (fT4l) is asymptotically stable. Let the input field have 
spectral density matrix R in [ico] . Then the output spectral density matrix is given by 



Rout[iuj} = E G [iu}R in [iu}E G [iuj^ . 



(41) 



If the input field state is vacuum p in = |0)(0|, then the steady-state output field state p out is 
a Gaussian state with covariance R ut(T~) which is in general not the vacuum state. However, if 
C+ = and Cl + = (passive systems) then we have 



R, 



out 



A(S G -H,0) 



/ 




A(S G -H,0) 



I 




(42) 



since for a passive system E G -[ico]E G -[iujY 
state is again the vacuum state. 



/ (Eq. (1241)). So in the passive case the output 
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III. Output intensities of quantum linear systems driven by multichannel 

PHOTON INPUTS 

We begin our study of the response of quantum linear systems to photon inputs by determining 
the statistics of the output field. Specifically, we consider multi-channel input signals, each 
channel with one photon, as defined in Section ITlI-Ai and then we find expressions for the output 
intensities (transient, Section Ull-B I and steady state, Section Ull-C I) and correlation, Section ITTl-DI 



A. Multichannel photon fields 

We now consider m field channels each of which is in a single photon state \ l Uk ), determined 
by possibly distinct pulse shapes satisfying the normalization condition \uk(t)\ 2 dt = 1, 
k = 1, . . . , m. This means that the state of the m-channel input is given by the tensor product 

i*„> = ii,,) ® • • • ® = n^K)r m >, (43) 

k=l 

where the m-channel vacuum state is denoted |0® m ) = |0) <g) • • • <g> |0). Here, Bl(i/k) = 
v k {t)b* k {t)dt is the creation operator for the /c-th field channel. 
For convenience, we let JFq denote the class of m-channel photon input field states: 

{™ POO ) 

\**) = l[B*k(yk)\0'* m ): / \v k (t)\ 2 dt = h * = U..,ffl • (44) 
fc=i ) 

B. Output intensity when the system is initialized at time to 

In this section, we study the intensity of output fields of the system G driven by the m-channel 
input field in the class (Eq. (|44)) ). We define this intensity to be 

n out (t) := (^ u \bL(t)bL(t)\^u) ■ (45) 



We first introduce some notations. For each k — 1, . . . , m, define 

|C fc ) := J] B*{v 3 W m ) = |U »---«> IWJ S |0) ® % k+l ) ® ■ ■ ■ ® \l Vm ). (46) 

Denote 

£- n {t) = dmg{v l {t),...,v m {t)), (47) 
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and 



C~ = [I m O m ]C, B~ = BS 



, B + = BS 



(48) 



For t > t , define classical variables 

(</>* v \tf(t)\<Ki) 



m_(t) := 



(0^|at(t)|0C m ) 



m+(t) :-- 



mx 2n 



(<K m \&(t)\<i>v v ) 



(49) 



mx In 



and 



£„(t) := (<py„\a(t)tf{t)\(py u )2nx2n, t > t . 



(50) 



We have the following result. 

Theorem 3: When the system G is driven by an m-channel photon input field \^ v ) in the 
class F Q (Eq. (144)) ). the output intensity n out (t) is given by 

noutit) =C~*^ u {t) T C~ T - C-*J n C~ T + St& n {tf&)S- (51) 
+ Si^m-{t)*C- T + C-#m„(t) T tr n (t)ST, 

where m_(t), m + (t), and E^i), defined in Eqs. (I49~l)-(l50l). satisfy the following differential 
equations: 

m_(*) = m-ft)^ + e^W 1 "^, t > t (52) 
m+(*) = m+(t)At + & n {t)B + \ t > t , (53) 

and 

± v (t) =AS v (t) + £„(t)At + fl-£;(t)m_(t) + B+C n (*) t m + (t) + m_(t)^r (t)tfi-t (54) 
+ m + (t)tf r (t)£+t + £?Sdiag(/ m , m ) < S , t jB t ) t > ^ 

respectively, with initial conditions m_(t ) = 0, m + (t ) = 0, and T, u (t ) = diag(/ m ,0 m ). 

Proof: It is straightforward to derive Eqs. (1521) and (1531) . To apply Lemma [2] to establish 
Eq. (|541) , it suffices to evaluate T(t) and F defined in Eqs. ( |36l) and (|37l) respectively. Notice 
that 

(<f>* u \a{t)tf{t)\<f>*„) = m_{t)^7 n (t)\ {^ v \a{t)b T {t)\^ v ) = m + (*) f C(*)- (55) 
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We have 



r(t)= m-W&W m + (tm n (t) 



On the other hand, for the single-photon input field, F defined in Eq. (1371) is 



Fdt = (y u \dB(t)d&(t)\V v ) 



Ira 





dt. 



(56) 



(57) 



Substitution of Eqs. (I56l)-(l57l) into Eq. (l38l) yields Eq. (l54l) . Finally, note that the gauge process 
of the output field A out (t) satisfies, e.g. US Sec. IV], 

dKutit) = S*dA(t)S^ + StdB*(t)L T (t) + L*(t)dB T (t)Sl + L*(t)L T (t)dt (58) 

= S*dA(t)££ + StdB*{t)a T {t)C- T + C-*a*(t)dB T {t)S T _ + C-*a*(t)a T {t)C- T dt. 

Noticing that c/A out = b* ut (t)b^ ut {t)dt and a(i)a l '(t) = (d(t)#d T (t)) T + J n , it can be readily 
shown that n out (t) defined in Eq. (l45l) satisfies Eq. (|5TI) . The proof is completed. ■ 

C. Output intensity in steady state 

In this subsection, we compute the steady-state output intensity when the system G is driven 
by an m-channel photon input field \9 V ) in the class F$ (Eq. (l44l)). 

The following is the main result of this subsection, whose proof is given in the Appendix. 

Theorem 4: The steady-state output intensity of the output fields of the system G driven by 
the m-channel single-photon input field \^f v ) in the class Fq (Eq. (l44l) ) is given by 



9G + (r)*g G+ (r) T dr + CM* Cut(t) T + Ct(*) # Ct(tf 



where 



/OO POO 
g G -(t-r)^ n (r)dr, g G +(t - r)& n (r)*dr. 

-OO J — OO 

In particular, the total output intensity is given by 

m m m 

Tr[rw(*)] = £ / \9^ + {t)\ 2 dt + £ |£f (t)| a + £ ICf*(*)| 2 , 



(59) 



(60) 



(61) 



j t k=l " u j,fc=l j,k=l 

where g G+ (t) is the element of g G + on the jth row and kth column. The same applies to g 3 ^{t), 

Sout ' dllu Sou* • 

For the single input case, the steady-state output intensity is given by 

\g G +(r)\ 2 dr + 



riout(t) 



9G+(t-r)C in (r)*dr 



9g- (t - r)£ in (r)dr 



(62) 
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Remark 4: For the single input case, if t = 0, in the Appendix we have shown that substi- 
tution of Eqs. (fTITT) . (fT22l) . (fl24l) and (fT27l) into Eq. (fT20l) yields 



/o 



Jo 



9G-(t-r)£ (r)dr 



(63) 



Assuming further that C + = f2 + = fL 
Eq. d63]) reduces to flU Eq. (26)]. 



0, C_ = «//^, and ^~(t) is v{t) in [|30l Eq. (18)], then 



D. Output covariance function in steady state 

If the input field is a multichannel photon field state \^> v ) as defined in Eq. (l43l . then it is 
easy to show that the input covariance function is 

R in (t, r) = ^ v [b(tfi(r)\ = R (t - r) + A(£(f), 0)A(£r( r ), 0)+, (64) 

where ^(t) is given in Eq. (|47T) and i? is the vacuum covariance defined by (|33~T) . According 
to Theorem [T] the output covariance function is 

/oo 
<7 G (t - r)di a g(/ m ,0 m )(fe(r - r)W + A(C t (t), ^(0)^(^(0, £ o + u t(0) f , 
-oo 

(65) 

with and £+ t given by Eq. Clearly 

-Rm(t, r) = i2in(r, t) 1 , -R ut(*, r) = R ou t(r, t) ] '. (66) 

Example 1. (Optical cavity) An optical cavity G is a single open oscillator lfl3l . 1131, ll40ll with 
Q_ = u £ R, Q + = 0, C_ = y/R, C + = JMl section IV. B.]). Let the pulse shape of the 
single-photon input field state | l v ) be given by 

x/2^e" 7 *, t > 0, 

" (67) 
0, t<0. 

The state |1„) can describe a single-photon field emitted from an optical cavity with damping 
rate y/Ty. The input covariance function is given by (l35l) . On the other hand by Eq. (l60l) . 

/2^ e -7* _ ( e -7* _ e -(f +&0*) , t > 



1/(0 



Ct(*) 



0. 



t < 



, Ct(*) = o. 



(68) 



Write 



_ ft . e -(f+^)(i-r) ; t > r 

5(t — r), t = r 

_ Ke -(f-Kw)(r-t) t<r 



(69) 
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It can be checked that the steady-state output correlation function is thus 

1 



Rout(t,r) = x(t,r) 







A(^(t),0)A(^(r),0)i 



(70) 



Clearly, the mean input intensity is n in (t) = \v(t)\ 2 . By Theorem HI the steady-state mean output 
intensity is n out (t) = |£~ ut (t)| 2 . According to Theorem [31 for t > 0, 



m„(t) 



2 - 7 - iu> 



-it _ e -(f-M* o 



m + (t) 



- 7 + icu 



Q g-7* _ e -(f -Hw)< 



(71) 



Hence, the evolution of system variables covariance £„(i) given in Eq. (l54l) can be expressed 
explicitly. In particular, when oj — 0, 

1 







27K 



(f-7) s 



( e -f* _ e" 7 ') 2 ^, t>0. 



(72) 



Example 2. (Degenerate parametric amplifier) A degenerate parametric amplifier (DPA) is an 
open oscillator that is able to produce squeezed output fields |[T3l . 0, ||4~0"1 . A model for a DPA 
is ttH pp. 220 and Chapter 10] 



a 



k — e 
— e k 

= v^a(t) + 6(t), (0 < e < «). 



(t)-V«£(f), a(0) = a, 



(73) 



If the single-photon input has the pulse shape defined in Eq. (l67l) . then by Eq. (l60l) . 



2 7 e- t7 + 



2K V / 27e~^ t ( e 2~~ *(27-k)-(27-k) cosh ^+esinh ^ 
(K-2 7 -e)(K-27+e) 



t > 
t < 



2K^/Tye ( -ee ^V^'+ecosh ^-(27-ft) sinh ^ 



(«-27-£)(ft-27+e) 



0. 



t > 
t < 



(74) 



(75) 



When the system is initialized at time t = 0, by Eq. (1631) . the mean output intensity is 

n out (t) = J* e— sinh 2 ^dr + \C at (t)\ 2 + |£+ t (*)| 2 , * > 0. (76) 
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Since |^ n ) = B*(v)\0), the covariance function corresponding to | &R m ) (the vacuum field 
state) is R (t) in Eq. (l33l) . According to Theorem [U the output covariance function is 

KR a {r) - Ke AT diag(l, 0) + diag(l, 0)5(r), 
KR a (r) - |/te AT diag(l, 0) - §«di&g(l, 0)e Atr + diag(l, O)tf(r) 
KR a {r) - Kdiag(l, 0)e" Atr + diag(l, 0)5(r), 



Rout,q\T~) 



T > 

r = , 
r < 



where 



e Ar T, 



T, 

Y e -At, 



r > 
r = 
r < 



T 



k / e 
'o 



At 



1 




(77) 



(78) 



As a result, according to Eq. (1651) . the covariance function of the output field of the DPA driven 
by the single-photon input field can be expressed as 



R out (t,r) = R^Jt -r) + A(Cut(t),Cut(t))^ t (r),Cut(r)) 



(79) 



where R ou t, g {t — r) is given in Eq. (1771) . 

Remark 5: In the cavity example, the output correlation function has the same form as the 
single photon correlation function (1351) . while this is not the case for the degenerate parametric 
amplifier. Thus the steady-state output state may in some cases be a single photon state, while 
in other cases more complex states may result. A general class of output states is the subject of 
the following section. 



IV. Photon-Gaussian states 

Since the states produced as the outputs of linear quantum systems need not necessarily be 
photon states, we consider a larger class of states that has the property that if a state in this 
class is input to a linear quantum system, then in steady state the output state is also in this 
class. This class of states is defined and studied in Section HV-BI However, the expressions for 
specifying these states are quite complicated, and so in Section IIV-AI we consider the simpler 
single channel case for pedagogical reasons. The calculations involved in determining the output 
states make use of the stable, but non-causal, inversions discussed in Section III-DI 
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A. The one channel case 

Given an initial joint system-field state p 0g = \<p) (<f>\ ® |0) (0|, denote 

Poog = lim U (t,t ) p 0g U (t,t )* . (80) 

t— >oo,to— >— oo 

Define 

Pfield,g = T^sys[Poog], (81) 

where the subscript "sys" indicates that the trace operation is with respect to the system. 
According to Theorem [2l pfi e id, g is the steady-state output field density with covariance function 
R out [icu] given in Eq. (I4TT) . 

We are in a position to prove a result concerning the output field of quantum linear systems 
driven by single-photon states. 

Proposition 2: Let m — 1 and suppose the input state p in = is a single photon state. 

Then the steady-state output field state for the linear quantum system G is given by 

Pout = (B*(Cut) ~ B(C llt ))Pfreld,g(B*(Cut) ~ £(&t))* (82) 

where 

A(CuM Cutis}) = Ea[8]A(u[s],0), (83) 



and pfieid, g , defined in Eq. (1811) . is the steady-state density operator for the output field with zero 
mean and covariance function 

Rout[M = E G [iuj}R in [iuj}E G [iuj] j , (84) 



where R in [iuj\ is the Fourier transform of Rq(t) defined in Eq. (1331) . 
Proof: The initial joint system-field density is 

p o = |0) (<t>\®\\ v )(l v \ = B*{v)p Qg B{v). 

The steady-state joint system-field density, denoted by p^, is 

Poo = lim U (t,t ) p U (t,to)* (85) 

t— >oo,to— >— oo 

= lim U(t,t )B*(v)p 0g B(is)U(t,t o y 

t— >oo,to— oo 

= lim U(t,t )B*(u)U(t,t y PoOtg lim U (t,t ) B(v)U (t,t )* , 

t— >oo,to— oo t— >oo,to— >— oo 
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t— >00, 



where p^g is given in Eq. < [8"0l ). Now we find expressions for the other terms. 

lim U(t,t )B*(v)U(t,t )* 

t— >OC,tQ— > — oo 

/oo 
v{r)I sys ®b*(r)drU{t,t o y 
-oo 

lim U(t,t ) [ v(r)I sys ®b*(r)drU(t,t )* 

O,t ->~CX3 

/OO 
u(r)b~*(r, -oo)dr 
-oo 

where 6" (Mo) = U(t,t )b(t)U*(t,t ). Now 

L 4 (t) = ^(t,to)^)^^o) = Ce A(t -* o) a+ / <7 G (*-r)6(r)dr 

•/to 

and so 

= Ce A{ *- to) a + / # G (t - r)b~(r, t )dr. 

J t 

Now send to — > — oo to obtain 

/t poo 
ga(t — r)b~(r, — oo) dr = / gait — r)b~ (r, —oo) dr. 
-oo J — oo 

Next, using the stable inverse of G (Lemma [Q Eq. (1261) ) this implies 

/oo 
g^- 1 it ~ r)b(r)dr 
■oo 



lim U(t,t )B*(v)U(t,t )* = I sys S (5*(C t ) - 

t— ^oo,£q— >— oo 



Therefore 



using Eq. (l26l) and the definitions of i out {t) and £^(2) in Eq. (|60l) . 
By Eqs. ([85]) and (gB, 

Poo = (/.„. ® - 5&)))poo,,(/^ ® (b*(c*) - 

This, together with Eq. (l8~TT) . 

Pout Tr^^^Jpoo] 

= Tr sys [(I sys ® (B*(C t ) - £(C + J))Poo, 9 (/^ ® - 

= (B*(Cut) - mD)pf*ei d AB*(Cut) - 

Eq. (f82l) is thus established. 
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The output state p out given by <[82l) is determined by a matrix £ oui = A(£~ ut , £^ t ) of functions 
which is obtained by convolving an input matrix £ in = A(v, 0) with the system G, and a Gaussian 
state Pfieid,g whose covariance R out is given by the usual transfer relation (|4TT) where i?j n is the 
covariance function for the vacuum field. It will be shown by Proposition [3] in Section [IV] that 
Pfieid, g in Eq. (f82l) is indeed normalized. 

Example 3. Refer to Example 1 on optical cavities. By Eq. ((42)) . pji e id, g is a vacuum state 
|0)(0|. By Proposition [2l the steady-state output field state is a pure state |\Pout) = B* (£~ ut )\0) , 
where is given in Eq. (f68|) . Clearly, the output is in a single-photon state. 

Example 4. Refer to Example 2 on degenerate parametric amplifiers. According to Proposition 
[21 the steady-state output state is p out = {B*{C ut ) ~ B(Cut))pfieid, g ( B *(Cut) ~ B(Cut))*, where 
iouti itut ^ given in Eqs. d74l)-(l75l). and the covariance function R out of Pfieid, g is given by Eq. 
(1771) . The normalization condition for p out will be given by Eq. (11001) . It can be verified that 
indeed Tr[p ou J = 1. Clearly, p out is not a single-photon state. 

B. Photon-Gaussian States 

In the previous section, in particular Example 4, we saw that a quantum linear system produces 
a somewhat complicated output state p out from a single photon input state. This output state was 
determined by pulse shapes and a Gaussian state. In this section we abstract the form of this 
output state and define a class T of photon-Gaussian states. The class J 7 contains single-photon 
states studied in ll29l and [|45l as special cases. Furthermore, in Theorem [5] we show that this 
class of states is invariant under the steady-state action of a linear quantum system, that is, 
p in G J= implies p out e J 7 . 

We first introduce some notations. Given t±, . . . , t 3 • G C and £i, . . . , Q G L 2 (C, C 2m ) where j 
is an arbitrary positive integer, define 



fi(-<i)®c---®c&(-*i)> 



ei(fi)®c---®c&(<j), 



^■(* 1 )®c-"®cCl(*i), 



(93) 



where ® c is the Kronecker product. Similarly, for the operators b(t), define 



M^t^j) :=b(t 1 )® c ---® c b(t j ). 



(94) 
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Finally for a matrix A, let A®? := A ® c • • • Cg> c A be a j-way Kronecker tensor product. Clearly, 
when j = 1, A®« = A. Define £(t) := A(f-(t), f +(t)) e C 2mx2w with entries £7 fc and £+ 0', = 
1, . . . , m) for matrices respectively. Let p R be a zero mean stationary Gaussian field 

state with correlation function R(r). 

The following equation will be used in Definition [Q 

/oo /*oo 
• • • / (M+(t^ m ) # ® c M 5 (t m+1 ^ 2m )) T J®" ® c e®"Tr[ j o Ji Mg(t 1 _, 2ro )]tZti . . . dt 2m = 1, 
oo */ — oo 

(95) 



2i/f 



where p# is a zero-mean Gaussian state with covariance function i?. 

Definition 1: A state p^R is said to be a photon-Gaussian state if it belongs to the set 

{mm /mm 
p<* = n E - Bit®.)) pr n E WW - 
fc=l j=l \fc=l j=l 

: £ and p fl satisfy Eq. (155} }. (96) 

Proposition 3: The photon-Gaussian states p^ R E T are normalized: Tr[p^ fi ] = 1. 
Proof: Partition to be 



r = ic 1 ■ ■ ■ rn e + = t^ 1 ■ • • e + ' m ], e = e 1 • ■ ■ e 

It is not hard to show that E T is in the form of 

Tr [P5,fi] 



2mn 



(97) 



(98) 



=Tr 



=Tr 



/OO / /"OO \ 

e fc (t) T e&(t)dt Pfl ] [ / e(tfob(t)dt 
\ I. J -oo I 



k • 

■ 

PR 



OO POO 



oo J — oo 



X 



/oo /*oo 
• • • / M S (t 1 ^ m ) T 6^Mg(t 1 _ >m ) C /t 1 • • • dt T 
■oo J —oo 



(M+(t^ m )* ® c M^t m+l ^ 2m )) T J^ ® c e®"Tr[ j o Ji M 8 (ti_, 2m )]tZti . . . dt 



2mi 



2 ni 
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OO POO 



[C(t)* -e(t)}Tr[p R b(t)b\r)} 



OO J — OO 



dtdr = 1. (100) 



where 6 is as introduced in the Notation part of the Introduction section. Thus, that Tr[p^ R ] = 1 
is equivalent to that Eq. (|951) holds. The proof is completed. ■ 
Note that when m = 1, a state p^R G T is 

ftlJl = (B*(C) - B(e)) PR (B*(C) - B(t)T, (99) 
and Eq. (1951) reduces to 

We are now ready to state the main result of this section. The proof is given in the Appendix. 
Theorem 5: Let p^ in ,R in G T be a photon-Gaussian input state. Then the linear quantum 
system G produces in steady state a photon-Gaussian output state p^ out ,R out G J 7 , where 

Ut[s] = E G [s]£in[s], (101) 
RouAM = E G [iu]R in [iuj]Z G [iuj] j! . (102) 

Without confusion, we may use the shorthand p out for P£ out ,R out . 

Remark 6: When the input state is a pure state and the system is passive (e.g., an optical 
cavity or a beamsplitter), it can be seen from Theorem \5\ that in steady state the output field is 
in a pure state. That is 

P^out,Rout I ^£out,Rout ) {^^out,Rout I ' 

In this case we use l^^,^^) to denote the steady-state output field state. Again, without 
confusion, we may use the shorthand \^ ou t) for \® z ou t,Rout) • 

We remark that the Gaussian part of the specification of photon-Gaussian states is needed to 
allow for quantum linear systems with active elements, such as degenerate parametric amplifiers, 
iTOl . Il40l . Il43l . In general, while passive devices produce vacuum from vacuum, active devices 
produce nontrivial Gaussian states from vacuum, and photon-Gaussian states from single photon 
states, as show in Example 3. 

This result provides a complete description of how quantum linear systems process highly 
non-classical photon-Gaussian states. In particular, the result provides the output response to 
single photon inputs. The result may be used for dynamical analysis, or for synthesis. Indeed, 
one could contemplate generalizations of the synthesis results given in [29] and in Section IVl to 
the class of photon-Gaussian states. 
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K) 




Fig. 3. A quantum linear stochastic system G driven by a single-photon state and a coherent state \a) simultaneously. 



Example 5. (Beam splitter) Beam splitters are archetype of static and passive quantum optical 
instruments 0, which can be modeled as b out (t) = S-b(t) where 

y/fj y/l-V 



s. 



with < 7] < 1. Let £r = diag(Vi, u 2 ), £^ = 0. According to Theorem \5\ and Remark [6l 

\Vout) 

pOO pOO 

= y/rj(l - rj) / ui(t)bl(t)dt v 2 {t)b\{t)dt\Qy) ® |0 2 ) 



(103) 



+7? / ^i(t)&J(t)dt|0i> g) / v 2 (t)b* 2 (t)dt\0 2 ) 

J —oo J — oo 

poo poo 

-(l-r/) / i/ 2 (t)6J(t)dt|0i) ® / i/i(t)6;(t)dt|0 



-^(l-^ld)® / VitybKQdt u 2 (t)b* 2 (t)dt\0 2 ). 



In particular, when 77 = 1/2 and = u 2 (t),\/t E R, the steady-state output state is 

v x {t)bl{t)dt ) |0 2 ). 



1 



l*out) = 2 



) ® |0 2 ) - ^|0i) ® 



(104) 



(105) 



In this case, the two photons can not exit from distinct output arms of the beam splitter. These 
results are consistent with the results of the calculations in flOl Sec. 16.4.2]. 

Example 6. (Linear quantum systems driven by both a single-photon state and a coherent state) 
The methods used in this paper may easily be adapted to treat multichannel input fields where 
some channels are coherent states while others are single photons. Consider a linear quantum 
system driven by a single-photon state \l v ) (channel 1) and a coherent state |a) (channel 2) 
simultaneously, Fig. [3] 
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Denote the pure input state by = \l u ) ® \a) = Bl(u) |0) £g> \a). It is easy to show that the 
correlation function for the composite state |0) <g) \a) is 

5{t-r) 

5(t-r) + a*(r)a(t) a{t)a{r) 



a*(t)a*(r) a*(t)a(r) 



R in (t,r) = (Oa|6(t)6+(r)|Oa) 



(106) 



Define 



(107) 



U(t)= gg-(t-r)v(r)dr, - r)i/(r)-dr, j = 1,2. 

By Theorem [5l it is not hard to show that the steady-state output state of the system driven by 
|^) is 



Pout = | - J p/w*« - ^(e^l)) 

where Pfieid, g has covariance function i? out given by Eq. (l40~l) . 



(108) 



V. Photon shape synthesis 

In this section we generalize the result of photon wavepacket shape synthesis results presented 
in ||29l . We consider a class of passive linear quantum systems for which C + = and Q + = 0. 
From Eq. (fl4l) it is easy to see that 

A ] + A + C ] C = 0, S ] C + (£5)+ = 0. (109) 

In this case E G [s] defined in Eq. (l22l) is £ G [s] = A(£ G -[s],0). By Proposition [T} 

~ G [-sf ~ G [s] = /, VsGC. (110) 

That is, H G [s] is all-pass [(471 pp. 357]. 

In what follows we study the following pulse shaping problem. Given two pulse shapes u(t) 
and u out {t) satisfying f™^ \v{t)\ 2 dt = f™^ \is out (t)\ 2 dt = 1, let u out [s] and i/[s] be the Laplace 
transform of v ou t(t) and v{t) respectively. If and u[s] satisfy 



\v\iuAr, Vcu e 



(HI) 



we show that, under mild conditions, there is an all-pass system which maps the input state \l u ) 
to produce the output field state \l Vout ). 
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Define 



!/[«] 



(112) 



We make the following assumptions on HgJs]: 
Assumptions 

1) is real-rational and has a state-space realization ll47l Chapter 3] 



A 


B 


C 


D 



(113) 



2) H(5 d [s] is Hurwitz stable. 

3) The above state space realization is minimal. 

4) D = I. 

Theorem 6: If the given input and output pulse shapes satisfy Eq. (11 1 1|) . and Ec d [s] satisfies 
Assumptions 1) - 4), then there exists a linear quantum stochastic system Eq[s] = A(Ea d [s],0) 
solving the pulse shaping problem: 



S_ = I, CL = C, C+ = 0, ft_ = - (XA - A^X) , n + = 0. 



(114) 



Proof: By Assumption 2, there exits a matrix X = X^ > such that 

A^X + XA + C^C = 0. (115) 

In fact, by Assumption 3, X > 0. Consequently, by Corollary 13.30 in [|47l . 

D t C + J B t X = 0, D j D = I. (116) 

Eqs. (II 151) and (II 161) can be rewritten as 

X" 1 ^ + AX" 1 + BB ] = 0, 5 = -X" 1 ^. (1 17) 

Then by Theorem 5. 1 in Il26ll . [s] can be implemented by an all-pass linear quantum stochastic 
system with parameters given in equation (II 141) . This completes the proof. ■ 
Example 7. Consider the following two functions v given by Eq. (1671) and 

2 (M^g-st _ 5(2 + t y-2t + i±p e (-i-^ ^ t > o 

0, t < 



Vout{t) 



(118) 
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Fig. 4. The upper plot is for v(t) and the lower plot is for i> ou t(i). The horizontal axes are time t, while the vertical axes are 
|^(£)| 2 (the upper) and |^ ou t(f)| 2 (the lower) respectively. 



It can be checked that v and v out satisfy Eq. (Ill II) . The shapes of u(i) and u out (t) are plotted 
in Fig. H Note that 

Vout[s\ = s-3 ^ s-(l-i) 

i/[s] s + 3 s + (l + i)' 

In the language of (£_, C_, C + , as discussed in Sec. III-AI by Theorem [6l the transfer 

function |=| corresponds to Gi = (l, -\/6, 0, 0, 0), while the transfer function corresponds 
G 2 = (l, -\/2, 0, 1, 0). Therefore, the whole system that transfers v to u out is a cascaded system 
made of G x and G 2 (TH Definition 5.3]. 

VI. Conclusions 

In this paper we have investigated the response of linear quantum systems driven by multi- 
channel photon input fields. Results concerning the intensity and correlations and states of output 
fields have been presented. In particular we have defined the class of photon-Gaussian states 
which arise when quantum optical systems with active components are driven by multi-channel 
photon input fields. Examples from quantum optics have been used to illustrate the results 
presented. Future work will include application of the results to specific problems in quantum 
technology. 
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Appendix 
Proofs 

Proof of Theorem S According to Eq. d20l) . 

n ut{t) = ( tt„, [ g G+ (t - r)*b(r)dr ( [ g G -(t - T)b(r)dr] * 

\ Jt \Jto / 

+ jf g G+ (t-r)*b(r)dr g G+ {t - r)b*(r)d?j m v 
+ U v , jf g G - (t - r)#b*{r)dr ( jf g G - (t - r)b{r)d^j tt„ 
j\ G -{t - r)*b*(r)dr (jf g G+ (t - r)6#(r)dr^ * 



(120) 



+ 



<0 W*0 

In what follows we evaluate each term of the right-hand side of Eq. (11201) . Firstly, it is easy to 
see that 



/ g G+ (t - r)*b(r)dT ( [ g G -(t - r)b(r)dr) 

k •/to \./to / 

¥„, jf <7g-(* " r)*b*(r)dr Qf g G+ (t - r)6#(r)dr) ^ 



0. 



0. 



Secondly, note that 



*»> ^ .9g+(*- 



(t - T) # &(r)dT ,g G+ (i - T)6 # (r)dT^ tf„ 



(121) 
(122) 

(123) 



m 



E£ij£*#(t-0**j(0& 

For given k and Z (k, I — 1, . . . , m), we have 



/ 9 G +(t - r)*g l ^+(t - r)dr + / 9 G +{t - r)*Uj(r)dr / g%+(t - r)vj{r)*dr. 

7=1 ^'o 7= 1 •'io ■'to 
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Consequently, Sending t — > — oo gives 

ft 



*„,J t g G+ (t-T)*b(r)dT^j g G+ (t-T)b*(r)dT^ tf„\ (125) 



^ G+ (r) # ^ + (r) r rfr + f+ t (t) # f+ t (*) 



Thirdly, note that 

ft 



/ g G -(t-r)#b*(T)dr( [ g G -(t-T)b(T)dr) V v 

J t a \Jt , 

Er=jto^-(*-o^w^ 



(126) 



EJLi J£ (* - Ofy ( L ) dL ■ ■ ■ E7=i II 9 G - (* - 6 i(0 d 



m pt ^mj 
3=1 J to 



For given k and / (k, I = 1, . . . , m), we have 



£j( flg.(*-r)*6J(r)drJ ^ jf ^_(f - rMOdrJ 

m /-t /-t 

/ g^{t — r)*i>j(r)*dr / g l ^_(t — r)vj(r)dr. 

7=1 "'*0 ■'to 



(127) 



As a result, sending £ — >■ — oo gives 



f 9 G -{t - T)#b#{r)dr ( fg G -{t - r)b(r)dr) ) = Ct(*) # 0*W T - ( 12 8) 

t ./to Vto / / 

Finally, substituting Eqs. (fTITT). (fl22l) . (fT23T) . and (fl28T) into Eq. (fl20l) yields Eq. @9J. Eq. ([6T) 
follows immediately from Eq. (1591) . This completes the proof. ■ 
Proof of Theorem^ First we prove that the steady-state output states are indeed in the form 
of Eq. (l96l) . In analog to Eq. (1851) . the steady-state joint system and output field state is 



lim tf(Mo)|$M®PUA»tf(Mo)' 



(129) 



^ ® II E J { {&n jk (r)bj*(r, -oo)dr - £^ fc (r)*&7 (r, -oo)) dr) Poo , £ 



k=i j=i 

m m 



X I 



sys 



II E y { (Cn Jfc (0^7*(r, -oo)dr - ^*(r)*67(r, -oo)) dr} 



fc=l 3=1 



where 



P 



oo,g 



lim f/(t,t o )|0)(0|® P/?in ?7(t,t o )*, 



t— >oo,trj — oo 



(130) 
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and for j = 1, ... , m, 

bJ(t,-oo)= lim U(t,to)bAt)U*(t,t ), bJ*(t,-oo)= lim U(t,t )b*Jt)U*(t,t ). (131) 
According to Eq. (l20l) . for any i e R, 

/■oo 

6(t) = / # G (t - r)U(r, t )b(r)U(r, t Q )*dr + e^-' )[/(t, t )aU(t, to)*. (132) 

Letting t — > — oo and substituting Eq. (11311) into Eq. (11321) yield &(£) = g G (t—r)b~ (r, —oo)dr. 
As a result 

/oo 
g G -i(t-r)b(r)dr. (133) 
-oo 

Partition g G - and g G + as 

<7g-(*) = [ g x G -{t) ■ ■ ■ g%-(t) ] , g G+ {t) = [ g G+ (t) ■ ■ ■ <#+(*) ] • (134) 
According to Lemma Q3 for j — 1, . . . , ra, 

/oo /*oo 
S G+ (i-r) t 6(0^+ / £ G _(t-r) T & # (t)^. (135) 
-oo J — oo 

In a similar way, 



6j (r, -oo) 



g J G -(L-r 



)%(i)dt 



g G+ (c-r) T b*(i)dc. 



By Eqs. (11351 ) and (11361) , for each j, k — 1, . . . , m, 



^n Jk (r)b;*(r, -oo)dr\<f><S> Rout ) - &*"*(r)*67(r, -oo)dr 



nn 



(136) 



(137) 



/ 9 3 G+ {i - r)%i)didr + &\t)I g^{i-r) T b#{i)didr 



m POO POO 

- E / st^w / » 

J —oo J — oo 



_{i - r)^b{i)didr + / Ct' W* / g j G+ {i-r) T b*(i)didr 



Substituting Eq. (11371) into Eq. (11291) gives 

^ ® n E(^(^ fe ) - ^(^f )) 
fc=i i=i 

Thus 



k=l j=l 



PUuR , t = ^sys [Poo] =I]E( 5 1^)- B i(^))^ ] IE( 5 ;(^f) " B Mti k )) 



k=l j=l 



Kk=l j=l 



(138) 
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where p Rout = Tr S yslPoo, g }, whose covariance function R out is given in Eq. (1102I ). Consequently, 
PHout,Rout is exactly in the form of Eq. ©. 

Next we prove that p^ ouU n out is normalized. Noticing that for a complex- valued function r}(i) 6 



r](t)b*(t)dt= / 77[-zo;]6*[iw]dw, / r]{t)*b(t)dt = / r/H*&Hdw. (139) 

) J -co J — oo J —oo 

It is easy to show that the frequency counterpart of p^ in ,R in in Eq. (1981 ) is 



Z71 



M ? (cji_ >m ) T e^="M g (wi_ >m )(ia;i^ m p/ ?i , 



x 



oo 



M+(w 1 _ >m ) t J^Mg(^ m )c/u; 1 _ 



where the shorthand dui^j is used to denote dw\ ■ ■ • duj for an arbitrary positive integer j. 
Consequently, similar to Eq. (1981 ), the normalization condition (1951 ) is equivalent to 



/oo /*oo 
■ ■ ■ / (M+(u^ m ) # ® c M ? -(w m+1 ^ 2m )) T J®™ ® c e®"Tr[p Rin M s (a; 1 ^ 2m )]da; 1 ^ 2m = 1. 
-oo J — oo 



2m 



(140) 



As a result, it suffices to show that p^ ouU R out has the normalization condition (11401) . Firstly, by 
Eq. (11011) . we have 



Partition £ out to be Ut = [Cut " ■ ' Then by Eq. 



^AP^outfiout] 



(141) 



2i/f 



xJ 8 ™ ® c e^TriPR^M^^du^. 
Secondly, noticing that ^ t [ico] = E G [iu]^ n [iu], Eq. (11411) becomes 

POO POO 

^[p^uRoJ = ■ ■ (M+Ju^)* ® c M- n ( Wm+1 ^ 2m )) T 



(142) 



2m 



x (S G [iw]t J)»? ® c (H G [-^] T 6)^Tr[p iJoiit M & (a; 1 ^ 2m )]^ 



► 2m,) 
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where slight abuse of notation is used, that is, 

(S G H* J)®" := (S G [f Wl ]t J) <g> c • • ■ ® c (S G [iw ro ]t J), (143) 

and 

(So[-iw] T e)^ := (H G [-za; m+1 ] T e) ® c • • • ® c (H G [-za; 2m ] T e). (144) 

Thirdly, denote 

S G H®= m := H G [zwi] ® c ■ • • ® c H G [zw 2m ]. (145) 

Noticing 

Tr^Mg^x^J] = ScH^^Tr^Mj^!^)], (146) 

Eq. (11421) becomes 

Tr[p fout ,* J (147) 



OO /"CO 

(M+ n (u;i^ m ) # ® c M- n (w m+ i_ >2m )) T 



OO J — OO 



2m 

x (H G H f JZ G [iu])®° ®c (S G [-za;] T 0H G H)®"Tr[p i?i . n M 6 (a; 1 ^ 2m )]da; 1 ^ 2m , 
where slight abuse of notation is used, that is, 

(H G H f JE G [iuj])®™ ■= (E G [iu^ JGiiui]) ® c • • ■ ® c (H G [zu; m ]VH G [iu; m ]), (148) 

and 

(2 G [-iw] T e2 G [iu;])®" 1 := (S G [-iw m+ i] T eE: G [?w m+1 ]) ® c • • • <g> c (H G [-iw 2m ] T e2 G [2u; 2m ]. 

(149) 

Finally, Eq. (11471) . together with the relations 

E G [iuj k ] j JE G [iuj k } = J, E G [-iuj k } T QE G [iuj k } = 6, (k = 1, . . . , 2m) (150) 
establishes Eq. (1140b . The proof is completed. ■ 
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